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The statistical properties of ensemble of disordered ID steric spin-chains (SSC) of various length 
are investigated. Using ID spin-glass type classical Hamiltonian, the recurrent trigonometrical 
equations for stationary points and corresponding conditions for the construction of stable ID 
SSCs are found. The ideal ensemble of spin-chains is analyzed and the latent interconnections 
between random angles and interaction constants for each set of three nearest-neighboring spins 
are found. It is analytically proved and by numerical calculation is shown that the interaction 
constant satisfies Levy's alpha-stable distribution law. Energy distribution in ensemble is calculated 
depending on different conditions of possible polarization of spin-chains. It is specifically shown that 
the dimensional effects in the form of set of local maximums in the energy distribution arise when 
the number of spin-chains M « iVj (where N x is number of spins in a chain) while in the case 
when M oc N% energy distribution has one global maximum and ensemble of spin-chains satisfies 
Birkhoff's ergodic theorem. Effective algorithm for parallel simulation of problem which includes 
calculation of different statistic parameters of ID SSCs ensemble is elaborated. 

PACS numbers: 71.45.-d, 75.10.Hk, 75.10.Nr, 81.5Kf 



I. FORMULATION OF THE PROBLEM 

Let us consider classical ensemble of disordered ID steric spin-chains (SSC), where it is supposed that interactions 
between spin-chains are absent (later it will be called an ideal ensemble) and that there are N x spins in an each chain. 
Despite some ideality of the model it can be interesting enough and rather convenient for investigation of a number 
of important and difficult applied problems of physics, chemistry, material science, biology, evolution, organization 
dynamics, hard-optimization, environmental and social structures, human logic systems, financial mathematics etc 
(see for example [3|— 18|] ) . As was shown by authors spin-glass model can be used for investigation of media's properties 
on scales of space-time periods of an external fields at conditions far from a usual equilibrium of media [9] . 

Mathematically mentioned type of ideal ensemble can be generated by ID Heisenberg spin-glass Hamiltonian 
without external field [IHl]: 

N x -1 

H (N X ) = - Ju+iSiS i+1 . (1) 

i=0 

where Si describes the i-th spin which is a unit length vector and has a random orientation. In the expression ([T]) 
Jii+i characterizes a random interaction constant between i and i + 1 spins, which can have positive and negative 
values as well fioj . 

In other words we consider the mathematical model of spin-chains ensemble where every spin-chain is like a regular 
ID lattice with the length L x — doN Xl where spins are put on nodes of lattice and interactions between them are 
random (see FIG 1). 

The distribution of spin-spin interaction constant W( J) is chosen from considerations of convenience and as a rule 
it is a Gauss-Edwards- Anderson model [Io| (see also [l|): 
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FIG. 1: A stable ID steric spin-chain with random interactions and the length of L x . The spherical angles <p and ip describe 
the spatial orientation of So spin, the pair of angles (ipi,ipi) correspondingly defines the spatial orientation of the spin Si, the 
distance between two neighboring spins in ID lattice is do- 



Let us recall that Jo and A J for this model are independent from the distance and scaled with the spin number 
N x as: 

(J) av = Jo cx N-\ AJociV- 1 / 2 , (3) 

in order to ensure a sensible thermodynamic limit. in Eqs. ([2]) and (j3]) describes the averaging procedure. 

Below we will investigate the issue of how much lawful the choice of this model is. 

For further investigations it is useful to rewrite the Hamiltonian (QJ in spherical coordinates (see FIG 1): 

Hq(N x ) — - ^2 Ja+i [cosipi cos-0i+i cos((fit - (fi+i) + sin -0; sin^+i] . (4) 
»=o 

A stationary point of the Hamiltonian is given by the system of trigonometrical equations: 

Oipi Oifi 

where 0j = {ij>i,ipi) are angles of i-th spin in the spherical coordinates system {tpi is a polar and ifi is an azimuthal 
angles), = (0i, 02-..-0jv a ,) respectively describe the angular part of a spin-chain configuration. 

Now using expression Q and equations ((3J) it is easy to find the following system of trigonometrical equations: 

i+l 

22 Jvi [ sm Vv - tan ipi cos tp v cos(ipi — <p v )] = 0, 

i/=i—l; y^i 
i+l 

Jvi cos ^ sin(</5j — ifv) = 0, J ui = J lu . (6) 

u=i—l\ v^i 

In case when all the interaction constants between i-th spin with its nearest-neighboring spins Ji-n, Ja+i and 
angle configurations (ipi—i, <Pi—i), i^iifi) are known, it is possible to explicitly calculate the pair of angles 0»+i = 
(i/>j+i,y>t+i)' Correspondingly, the i-th spin will be in the ground state (in the state of minimum energy) if in the 
stationary point 0° = (ip®, ip®j the following conditions are satisfied: 

^h* t (®°) > 0, A^ t (0°) A ViH3i (0°) - (0°) > 0, (7) 
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where A QiQi (9°) = d 2 H /da 2 , A aiPi (Q9) = A fta( (6?) = d 2 H /da t d(3 l , in addition: 



i+l . 

22 J»i [cos -01/ cos(<£>„ — + tan-0 4 ° sin0v] > cos^ , 

v=i— 1; i/^i ^ 



i+l 



Aw(®i) = ^ /J ^i/i COS "0!/ COS^,, -Ifi) >C0S^° 
u—i — 1 ; u^i 

A*l>t<t>i = ] X! ^ cos^r/ sin(^ - <£°) > sin-0 4 c 



(8) 



y — ?' — 1 ; v^i 

Taking into account the second equation in © we can reduce condition ([7]) to the following kind: 

^(e°)>o, A Vi ^(6?)>0. (9) 

So, with the help of Eq.s (|6]) and conditions ([9|) huge number of stable II? SSCs may be calculated and on its basis 
it is possible to further construct the statistical properties of ID SSCs ensemble. It is important to note that the 
average polarization of ID SSCs ensemble is supposed to be equal to zero. 

Now we can construct the distribution function of energy in ID SCCs ensemble. To this effect it is useful to divide 
the nondimensional energy axis e = e/de into regions > sq > ... > e„, where n >> 1 and e is a real energy axis. The 
number of stable ID SSC configurations with length of L x in the range of energy [e — 5e,e + 5e] will be denoted by 
Mi x (e) while the number of all stable ID SSC configurations - correspondingly by symbol M[ ul1 — ^L x {^j)- 
Accordingly, the energy distribution function into the ID SSCs ensemble may be defined by expressions: 

F Lx (s;d (T))=M Lx (e)/M[f, (10) 

where distribution function is normalized to unit: 

n „o 
hm y2F Lx ( Sj ;d (T))S Sj = / F L „(e;d (T))de = 1. 

By similar way we can define also distributions for polarization and for a spin-spin interaction constant. 

II. ALGORITHM OF ID SSCS IDEAL ENSEMBLE SIMULATION 

Now our aim is elaboration of algorithm for parallel simulation of ideal ensemble of II? SSCs. 
Using equations ([6]) for stationary points of Hamiltonian Hq(N x ) we can find the following equations system: 

Ji-iifsinV'i-i - tanV'i cos V'i-i cos(<pi - ifi-i)] + Ji i+i [sin 0>j + i 

- tan0>i cos0 i+ i cos(<^ - tfi+i)] = 0, 
Ji-n cos0> 4 _i sm(ipi - tpi-i) + Ji+n cos?Ai+i sin(^ - <p i+ i) — 0. (11) 

After designations: 

x = cosijji+i, y = sm(ipi-(p i+ x), (12) 
the system (|llj) may be transformed to the following form: 



Ci + J ii+1 [y/l - x 2 - tam/^aVl - y 2 ] = 0, C 2 + Ja+i xy = 0, (13) 
where parameters C\ and Ci are defined by expressions: 

Ci = Ji-ulsmipi-! - tan^cosV'i-i cos(( ( c i - ipi-i)], 

C 2 = Ji-ncosipi-i sin(cpi - ifi-x). (14) 
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From the system of equations (|13[) we can find the equation for the unknown variable y: 



Ciy + CaVT^tan^ + ^J 2 l+1 y 2 - Cf = 0. (15) 
We can transform the equation (p~5|) to the following equation of fourth order: 

[A 2 + 4CfCi sin ^) y 4 - 2 [AC\ + 2dCf sin 2 V;] ?/ + C\ = 0, (16) 

where 

A = Jfi+ 1 cos 2 i>i - C\ + C 2 2 sin 2 ^ . (17) 
Discriminant of equation (|16|) is equal to: 

D = C%(A + 2Ci sin 2 ^) 2 - C 2 4 (v4 2 + 4C 2 C* 2 sin 2 ^) 

= 4C|C 2 sin 2 ipi ( A + C\ sin 2 ip, - C|). 

From the condition of nonncgativity of discriminant D > we can find the following condition: 

A + C\ sin 2 Vi — C| > 0. (18) 

Further substituting the value of A from (JTTJ) into (|18[) we can find the new condition to which the interaction constant 
between two successive spins should satisfy: 

^ +1 >C 2 + Cf. (19) 
Now we can write the following expressions for unknown variables x and y: 

x 2 ° 2 



J2 2 ' 

cos 2 in J 2 l+1 ± 2Ci sin ipi cos ipiJ Jf i+1 - C\ - C 2 + C 3 + 2C 2 sin 2 ^ 

v 2 = C 2 - (20) 

2 cos 4 ^j4 +1 +2C 3 cos 2 ^J 2 l+1 + (C 1 2 +sin 2 V l C 2 ) 2 ' V ' 

where C 3 = -C\ + C\ sin 2 ^. 

Finally taking into account designations (|12j) we can find new conditions of restriction of the calculated angles 
((p i+ i,ip i+ i): 

< x 2 < 1, < y 2 < 1. (21) 
These conditions are very important for elaborating correct and effective algorithm for numerical simulations. 

A. Algorithm description 

This is parallel algorithm for simulation of ID SSCs ensemble, which consists of separate iterative calculations of 
nodes in ID SSC. The first and second nodes are initialized randomly, then i-th node is obtained from (i — 2)-th and 
(i — l)-th layers nodes. Every node contains the following information: 
cp-polar angle, 
-0-azimuthal angle, 
J-interaction coefficient, 

The following parameters are initializes in the following way: 
ipo and <fii - rand()*2*7r*i?; 
ipo and ipi - acos (rand()); 
J i - rand(); 
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where rand() function generates uniformly distributed random numbers on the interval (0, 1). 

The algorithm pseudo-code is following: 

// generate n separate independent sets of problem in parallel 
for i = 1 : N x 

for j = 1: R If regenerate Ji maximum R times if needed 

for k = 1 : Li //go through all elements in the i-th layer if conditions 
// © a re satisfied 

begin 

/ / calculate energy on i-th layer, 

/ / calculate polarization on x, y and z-axis 

// calculate and 

/ / save Ji value 

end 

cndfor 

endfor 
endfor 

if (i == N x ) If reached the iV x -th layer 
begin 

/ / save energy, polarizations values 

end 

cndif 

// construct distribution functions of energy e, polarization p and 
// interaction constant J 

II calculate the mean value of energy e, polarization p, interaction constant J and 
// its variance J 2 . 

III. NUMERICAL SIMULATION 

We will consider an ideal ensemble of ID SSCs which consists of M number of spin-chains each of them with the 
length 25do- For realization of parallel simulation we will use algorithm A (see FIG 2). 

The parallel algorithm works in the following way. Randomly M sets of initial parameters are generated and 
parallel calculations of equations (|20[) for unknown variables x and y transact with taking into account conditions 
(1211) . However only specifying of initial conditions is not enough for solution of these equations. Evidently these 
equations can be solved after definition of the constant Jqi, which is also randomly generated. In the case when 
solutions are found then conditions of stability of spin in node © are checked. The solution proceeds for the following 
spin if the specified conditions © are satisfied. If conditions are not satisfied, a new constant Jqi is randomly 
generated and correspondingly new solutions are found which are checked later on conditions ©• This cycle on each 
spin repeats until the solutions do not satisfy to conditions of the minimum spin energy in the node. 

At first we have conducted numerical simulation for definition of different statistical parameters of the ensemble 
which consists of 10 2 spin-chains. Let us recall that the number of simulation of spin-chains define the number of 
spin-chains in the ensemble. As the simulation shows (see the left picture in FIG 3) the energy distribution function 
has a set of local maximums (e (0) ,..., sM). Obviously they are dimensional effects and are similar to the first-order 
phase transitions which often happen in spin- glass systems 

Let us note that during simulation we suppose that spin-chains can be polarized up to 20 percent i.e. the total 
value of spins sum in each chain can be in an interval of — 5 < p < 5, where p designates the polarization of 
spin-chain. In other words each spin-chain is a vector of certain length which is directed to coordinate x. As 
calculations show, in the ensemble consisting of a small number of spin-chains, for example, of the order 10 2 , the 
self-averaging of spin-chains does not occur in full measure i.e. the total polarization of an ensemble differs from zero: 
p x = —0.33099, p y = —0.035191, p z — —0.024543 where p = f_ °° F{p)dp, where it is supposed that p — (p x , p y , p z ). 

In this case the average energy of an ensemble is equal to e = —14.121, where e = f_ F(e)ede. 

For the ensemble which consists of 2.10 3 spin-chains (see FIG 4), the dimensional effects practically disappear. The 
summary polarization of ensemble in this case is very small: p x = —0.020538, p y — —0.047634, p z — —0.12687 and 
correspondingly the average energy of ID SSC is equal to e = —13.603. 
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FIG. 2: The algorithm of ID SSCs of ideal ensemble parallel simul^^on^^^^^^^^^^^giieters. 
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FIG. 3: The energy distribution where there are apparently many local minimum of energy for ensemble of ID SSCs with the 
length of L x = 25<io, which consists of 10 2 spin-chains (the left picture). On the right picture polarization distributions of 
ensemble on coordinates x,y and z are shown. 



Ensemble which consists of 10 4 spin-chains has an energy distribution F(e) with one global maximum (see Fig 5). 
As to polarization distributions, F(p x ) F(jp y ), and F(jp z ), in the considered case are obviously very symmetric in com- 
parison with similar distributions of previous ensembles (see FIG 3 and Fig 4) . The average values of polarizations on 
coordinates for this ensemble are much smaller p x = —0.0072863, p y = —0.014242, p z = —0.018387, correspondingly 
the average energy is equal to e = —13.634. Thus in the case when ensemble consists of a big number of spin-chains, 
the self-averaging of spin-chains system occurs with high accuracy. Whereas the summation procedure on the number 
of spins in chain or spin-chains ensemble is similar to the procedure of averaging by the natural parameter or " timing" 
in the dynamical system, it is possible to introduce the concept of ergodicity for the both separate spin-chains and 
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FIG. 4: In the left picture is shown the energy distribution in the ensemble of ID SSCs with the length of L x = 25do, which 
consists of 2-10 3 spin-chains. Apparently, the number of local minimums of energy is promptly reduced comparing with the 
increase of spin-chains. On the right picture polarization distributions of ensemble on coordinates x, y and z are shown. 




FIG. 5: The energy distribution and its fitted curve (left picture) in ensemble of ID SSCs with the length of L x — 25do, 
which consists of 10 4 spin-chains. Evidently there is only one global maximum for energy distribution. In the right picture 
polarization distributions are shown correspondingly on coordinates x,y and z. 



ensemble as a whole. 

Thus as calculations show Birkhoff ergodic hypothesis [ll[ may be used for ensembles which consist of M ~ N% spin- 
chains in order to change the summation of spin-chains on the integration by the energy distribution of the ensemble. 
The energy distribution of ensemble does not depend on the length of the spin-chain in the limit of ergodicity and it 
can be fitted very precisely with Eckart function [l2[ (see FIG 5, the smooth 

F(e) = C(o, b, c, 7 ) ( - „ + ^- 1 , (22) 




FIG. 6: The energy distributions for ensembles consisting of ID SSCs of the length L x = 25do, with spin-chains polarization 
correspondingly up to 20, 40 and 100 percents (left picture). Note that all the ensembles consist of 10 4 spin-chains and their 
distributions practically do not differ. On the right picture the distribution of the spin-spin interaction constant is shown which 
differs essentially from Gauss-Edwards- Anderson distribution model (2). 



where a, b, c and 7 some constants, in addition C is a normalization constant and can be found from the condition: 

i-o 



I F{e)de = 1. (23) 

J — OO 



C- 1 (a,6,c, 7 ) = ^ln(l + &) + ^. (24) 



By placing ([22j) into ([23)) we can find: 



After fitting the energy distribution by means of analytical function (1221) we find values of constants by entering into 
the function: a = 131.4, b = 3138.2, c = -1.20344 and 7 = 0.162174. 

We have also calculated ID SSCs ensemble with the length of spin-chains 25g?o and correspondingly with polariza- 
tions of spin-chains up to 20, 40, and 100 percents (see Fig 6, the left picture). In particular, as it follows from the 
picture the energy distribution does not depend on the degree of spin-chains polarization. We also have conducted 
simulation of ensembles which consist of spin-chains with lengths 100<io and lOOOdo correspondingly. As the numer- 
ical modeling shows, statistical properties of ensembles are similar. In the considered cases distributions of energy 
concentrate correspondingly on scales 100c£o and lOOOdo- Limits of ergodicities of ensembles are also investigated and 
it is shown that in these cases too it is of an order N%. 

Finally it is important to note that the distribution of spin-spin interaction constant is not defined apriori with 
the help of expression ^ but with the mass calculations of equations ([5]). On the basis of the obtained numerical 
data, the distribution of interaction constant W(J) = F(J) is constructed (see Fig 6, the right picture) from which it 
follows, that it essentially differs from the Gauss-Edwards- Anderson distribution model ((21). The obtained distribution 
relatively is well fitted by the normalized to the unit of nonsymmetric Cauchy function [13J: 

n[g 2 + {J- aoy\ 

where g, /3 and do are some adjusting parameters which are found from the condition of a good approximation of 
the data visualization curve. In the considered case they are correspondingly equal to: g = 0.27862, f3 = 0.009 and 
ao = 0.083236. Nevertheless, as the detailed analysis of curve of numerical data visualization shows (in particular 
its asymptotes) the distribution of interaction constant can be approximated precisely by Levy skew alpha-stable 
distribution function. Let us recall that Levy skew alpha-stable distribution is a continuous probability and a limit 
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of certain random process X(a,/3, 7, S;k) where parameters describe correspondingly: an index of stability or char- 
acteristic exponent a £ (0;2], a skewness parameter (3 6 [—1; 1], a scale parameter 7 > 0, a location parameter S € K 
and an integer k shows the certain parametrization (see in more detailed references [H], El). Let us note, that the 
mean of distribution and its variance are infinite. However, taking into account that spin-spin interaction constant 
has limited value in real physical systems, it is possible to calculate distribution mean and its variance. In particular 
if J 6 [-5, +5] then 7 = 0.50113 and 7 1 = 2, 1052. 

IV. CONCLUSION 

The investigation of statistical properties of classical spin-glass system of various sizes is very important for un- 
derstanding possibilities of effective influence and control over parameters of medium with the help of weak external 
fields. Evidently, when we put the spin-glass in external field the space-time periods define scales on which prob- 
ably an essential changes in medium occur. For simplicity we suppose that the spin-glass system is an ensemble 
which consists of disordered ID steric spin-chains of L x lengths, between which interaction is absent (ideal ensem- 
ble). This type of classical ensemble is described by Heisenberg Hamiltonian ©. We have researched conditions of 
arising of stable spin-chains Eqs. (fTTj) and nonequalities © and found a latent connection between random variables 
(see expression (|19l) ). which shows that the distribution for spin-spin interaction constant can not be described by 
Guss-Edwards- Anderson model. In the result of equations of stationary points analysis ([lip we have found system 
of recurrent equations (|2"0"]l and new conditions (|2"Tj) . On the basis of obtained mathematical formulas the effective 
parallel algorithm for numerical simulation is developed which was realized on the example of the ensemble which 
consists of ID SSCs with length 25<io- Similar to the dynamical systems, we have introduced the idea of Birkhoff 
ergodic hypothesis [ll[ for the statical spin-glass systems. In this case the number of spin-chains of ensemble plays a 
role of the natural or "timing" parameter of the system. Numerical simulations show that the ergodic hypothesis may 
be used for the case when ensemble consists of M oc N% spin-chains in order to change the summation of spin-chains 
on the integration by the energy (polarization, etc.) distribution of the ensemble. 

In particular, we have made numerical experiments for ensembles which include 10 2 , 2 • 10 3 and 10 4 spin-chains. 
As it was shown by simulations in the case when M <C N% for an ensemble, they are characteristic dimensional 
effects in energy distribution (the left picture on FIG 3). When the number of spin-chains is M of order 2 • 10 3 or 
more 10 , dimensional effects disappear and correspondingly energy distribution functions have one global maximum 
(see left pictures on FIG 4 and FIG 5). As it was shown, when increasing spin-chains number, the total and partial 
polarizations of the ensemble disappear. Let us note, that at modelling by algorithm (see scheme on FIG 2) condition 
(1191) specifies the region of localization of random interaction constant Ja+\ which depends on angular configurations 
(i — l)-th and i-th spins and interaction constant Ji-n between them. As a result, it allows to accelerate calculations 
of each spin-chain and hence the speed of parallel calculations of ensemble is increased essentially. 

Finally it is important to note that it is proved, that the spin-spin interaction constant Ja+i has a form of Levy 
skew alpha-stable distribution (see the right picture on FIG 6). The considered scheme of solution of ID steric spin- 
glass problem can be used in different applied fields (see e.g. It can also be useful for analyzing 3D spin-glass 
problem and creation of an effective parallel simulation algorithm of the spin-glass system with large dimensionality. 
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